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摘要 :【 目 的] 草地 贪 夜 蛾 Spodoptera frugiperda 是 一 种 新 近 入 侵 我 国 的 重要 害虫 。 本 研究 旨 在 对 草 
地 贪 夜 蛾 雄性 成 让 和 5 龄 幼虫 两 个 不 同 发 育 阶 段 的 转录 组 进行 比较 分 析 。【 方 法 】 利 用 高 通 量 测 
序 技术 对 草地 贪 夜 蛾 雄性 成 下 和 5 龄 幼虫 进行 转录 组 测序 和 数据 组 装 ,并 对 转录 组 数据 进行 功能 
注释 和 比较 分 析 。 【结果 】 经 de novo 组 装 共 获得 209 002 条 转录 本 ,平均 长 度 为 687.55 bp ,N50 为 
982 bp。 共 有 46 198 4& (57. 43% ) unigene 在 至 少 一 个 数据 库 中 获得 功能 注释 ,其 中 1713 条 
(2. 13% ) unigene 在 所 有 数据 库 中 均 能 获得 注释 。 在 GO 数据 库 中 获得 205 269 条 unigene 的 注 
释 ,主要 包括 68 个 功能 分 类 ;在 KEGG 数据 库 中 共有 3 408 条 unigene 得 到 注释 ,涉及 277 个 代谢 
通路 。 共 鉴定 到 424 个 噢 觉 相 关 的 基因 ,并 且 在 雄性 成 忠和 5 龄 幼虫 之 间 的 表达 存在 差异 。 通 过 
比较 转录 组 分 析 , 在 雄性 成 虫 中 鉴定 到 9 162 个 上 调和 6 399 个 下 调 差 异 表达 基因 (DEGs) ;功能 富 
集 分 析 发 现在 上 调 DEGs 中 涉及 信息 素 以 及 信号 转 导 的 代谢 通路 显著 富 集 , 而 下 调 DEGs 中 涉及 解 
毒 相 关 的 通路 显著 富 集 。【 结 论 ] 这 些 转录 组 数据 为 探究 革 地 贪 夜 蛾 的 生长 发 育 、 嗅 觉 相关 功能 基 
因 以 及 候选 分 子 靶 标 提供 资源 信息 。 

关键 词 : 草地 贪 夜 蛾 ; 发 育 阶 段 ; 嗅觉 基因 ; 转录 组 ; 基因 注释 ; 生物 信息 学 
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Abstract: [ Aim] The fall armyworm, Spodoptera frugiperda, is a newly invaded important pest in 


China. This study aims to comparatively analyze the transcriptomes of two developmental stages, male 
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adult and the 5th instar larva, of S. frugiperda. [Methods] High-throughput sequencing technology was 
used to sequence and assembly the transcriptomes of male adult and the 5th instar larva of S. frugiperda, 
and the functional annotation and comparative analysis of transcriptome data were performed. [ Results] 
After a de novo assembly , a total of 209 002 transcripts were obtained, with an average length of 687.55 
bp and an N50 of 982 bp. A total of 46 198 (57. 4396) unigenes were annotated in at least one 
database, among which 1 713 (2.13% ) unigenes were annotated in all databases. In the GO database, 
205 269 unigene annotations were obtained, mainly including 68 functional classifications. A total of 
3 408 unigenes were annotated in KEGG database, involving 277 metabolic pathways. A total of 424 
olfaction-related genes were identified and showed differential expression in male adults and the 5th instar 
larvae of S. frugiperda. By comparative transcriptome analysis, 9 162 up-regulated differentially 
expressed genes ( DEGs) and 6 399 down-regulated DEGs were identified in male adult. Functional 
enrichment analysis revealed that the up-regulated DEGs were significantly enriched in the pathways 
involved in pheromone and signal transduction, while the down-regulated DEGs were significantly 
enriched in the pathways related to detoxification. [ Conclusion] These transcriptome data provide some 


valuable information for exploring the functional genes related to the growth and development and olfaction 


and candidate molecular targets in S. frugiperda. 
Key words: Spodoptera frugiperda; 


annotation; bioinformatics 


草地 贪 夜 蛾 Spodoptera furgiperda , X WTgoRs Œ, , 
Ja REH H (Lepidoptera ) T EE] ( Noctuidae) , 原 产 于 
美洲 热带 和 亚热带 地 区 , RAEE DOS) ^ XE KE 
强 食性 杂 ,为 害 重 和 防 探 难度 大 等 特点 ,是 广泛 分 
布 于 美洲 各 国 的 重大 迁 飞 性 害虫 (Sparks，1980 ) 。 
在 美洲 ,草地 贪 夜 蛾 已 分 化 为 玉米 型 和 水 稻 型 ,而 入 
侵 中 国 的 草地 贪 夜 蛾 属 玉米 型 ( 张 夭 等 ,2019 ) 。 
2019 年 1 月 ,在 我 国 云南 普洱 和 德 宏 市 首次 发 现 草 
地 贪 夜 峨 入 侵 (https: // www. ippe. int/en/news/first- 
detection-of-fall-armyworm-in-china/)。 草 地 贪 夜 蛾 
在 我 国 为 害 的 作物 有 玉米 .上 甘 苔 .上 甘蓝、 马铃薯 和 营 
改 等 ,但 玉米 比 其 作物 更 容易 受到 攻击 ,其 为 害 率 严 
的 可 达 80% ( 雇 永 林 等 , 2019)。 截 止 2019 年 7 
月 ,该 虫 已 扩散 至 全 国 20 个 省 (自治 区 )1 128 个 县 
(市 .区 ) ,累计 发 生 面积 约 831 71 Bi ( http: / www. 
sohu. com/a/325385834 114960) 。 

近年 来 迅速 发 展 的 高 通 量 测序 技术 可 为 昆虫 分 
子 生物 学 研究 提供 基因 数据 。 昆 虫 的 生物 学 性 状 表 
达 和 行为 密切 关联 ,基因 信息 可 以 有 效 指 导 筛 查 发 
现 和 研究 与 目的 性 状 表 达 、 遗 传 调控 相关 的 关键 基 
以 及 预测 调控 蛋白 的 功能 ,筛选 得 到 可 应 用 于 害 
Th EP BÉ EU dee D] ,为 草地 贪 夜 蛾 等 农业 害虫 的 防 
治 提供 新 的 思路 。Kakumani 55 (2014) 利用 第 二 代 
测序 技术 对 草地 贪 夜 蛾 卵 梨 分 离 细胞 系 (SF21) 进 
行 了 基因 组 测序 ,组 装 了 37 243 个 scaffolds, N50 为 









































developmental stage; 
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53.7 kb ,基因 组 总 大 小 为 3538 Mb ,其 中 350 Mb 的 基 
因 组 序列 比 对 到 了 28 RRE Bombyx mori 的 染色 体 
上 。2017 年 ,科学 家 们 还 分 别 对 玉米 型 和 水 稻 型 草 
地 贪 夜 蛾 进行 了 重 测序 ,发 现 了 一 些 基因 家 族 在 两 
个 品系 间 存 在 差异 扩 增 ,可 能 与 适应 不 同 寄主 植物 
的 进化 相关 ( Gouin et al., 2017)。Liu 等 (2019 ) 对 
入侵 中 国 的 草地 贪 夜 蛾 个 体 进行 了 全 基因 组 测序 ， 
结果 表明 中 国 的 草地 贪 夜 蛾 包含 玉米 型 和 水 稻 型 两 
个 品系 ,最 有 可 能 来 源 于 非洲 。 本 研究 采用 高 通 量 
测序 技术 ,初步 建立 草地 贪 夜 蛾 雄性 成 虫 和 5 龄 幼 
虫 发 育 阶段 转录 组 数据 库 , 为 解析 草地 贪 夜 蛾 不 同 
发 育 阶段 生物 学 信息 .挖掘 生长 发 育 关键 的 调控 基 
因 黄 定 基础 。 


















































1 材料 与 方法 
1.1 iX 
草地 贪 夜 蛾 在 广东 省 农业 科学 院 植物 保护 研究 
所 内 饲养 ,于 2019 年 6 月 份 取样 ,选取 健康 雄性 成 
THURIS 龄 幼虫 各 3 K, TENDRE SET SEA X 
氮 处 理 10 min ,随后 置 于 超 - 80% 冰箱 保存 。 
1.2 总 RNA 提取 

草地 贪 夜 蛾 总 RNA 采用 Total RNA Extractor 
(B511311-0025 ) 试剂 盒 ( 上 海 生 工 生物 工程 股份 有 
限 公 司 ) 进行 提取 ,然后 通过 1% 电泳 凝 胶 检 测 提取 
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的 RNA 的 完整 性 。 采 用 Invitrogen Qubit@ 2.0 X 
光 试 剂 盒 (Fluorometer，Life Technologies Invitrogen , 
Q32886) 对 总 RNA 进行 定量 。 
1.3 ”转录 组 测序 与 拼接 组 装 

委托 上 海 生 工 采 用 Hllumina HiSeq2500 高 通 量 
测序 平台 进行 转录 组 测序 。 测 序 得 到 的 原始 数据 ， 
通过 NGS QC 进行 质量 评估 和 Trimmomatic 进行 质 
量 剪 切 ,过 滤 掉 接头 、 低 质量 的 序列 (长 度 小 于 35 nt 
的 reads) ji N 碱 基 的 序列 和 低 质 碱 基 (Q 值 <20 ) 
的 序列 ,得 到 高 质量 的 clean data ( Patel and Jain, 
2019)。 使 用 Trinity 软件 对 clean data 进行 de novo 
拼接 组 装 ,再 采用 RSeQC 软件 去 除 转录 本 中 的 匈 余 
序列 ,得 到 unigenes( Haas et al., 2013) 。 
1.4 基因 功能 注释 

采用 blast 将 组 装 的 unigene 与 CDD, KOG, COG, 
NR, NT, PFAM, Swiss-Prot 和 TrEMBL 等 多 个 数据 
库 比 对 得 到 功能 注释 信息 ,根据 转录 本 与 数据 库 
blast 比 对 结果 和 transdecoder 进行 CDS 预测 。 根 据 
unigene 与 Swiss-Prot 和 TrEMBL 的 注释 结果 得 到 
GO 功能 注释 信息 ,利用 KAAS 得 到 KEGG 注释 信息 。 
1.5 基因 表达 分 析 

使 用 Salmon 对 unigene 的 表达 量 水 平 进行 估 
计 , 利 用 每 百 万 条 reads 的 转录 本 ( transcripts per 
million, TPM) 计量 在 RNA 池 中 某 个 转录 本 的 比例 。 
使 用 DESeq2 进行 基因 表达 差异 分 析 , 筛选 差异 表 
达 基 因 的 标准 为 FDR 小 于 0. 05 且 差 异 倍数 (fold 
change, FC) 大 于 2, 使 用 topGO 进行 GO 富 集 分 




































































析 ,使 用 clusterProfiler 进行 KEGG 通路 和 KOG 分 类 
富 集 分 析 。 








2 结果 


2.1 转录 组 测序 与 de novo 组 装 

草地 贪 夜 蛾 雄性 成 虫 和 5 龄 幼虫 cDNA. 文库 的 
构建 由 上 海牛 工 完 成 ,并 利用 Ilumina HiSeq 2500 
测序 平台 测序 ,雄性 成 虫 共 获得 46 851 886 条 raw 
reads , 总 碱 基数 为 7 027 782 900 bp; 5 龄 幼虫 共 获 得 
48 160 250 条 raw reads ,总 碱 基数 为 7 224 037 500 
bp。 使 用 Trimmomatie 对 原始 测序 数据 进行 处 理 ， 
去 掉 含 有 带 接 头 和 低 质量 的 序列 ,雄性 成 虫 共 得 到 
44 922 812 clean reads ,总 碱 基数 目 为 6 444 354 229 
bp, GC 含量 为 44.12% , Q30 值 达到 94.01% ;5 龄 幼 
虫 共 得 到 46 414 162 clean reads, 总 碱 基 数目 为 
6 679 121 985 bp,GC 含量 为 46. 40% , Q30 值 达到 
94. 5396 ,表明 文库 构建 质量 良好 ,测序 得 到 的 数据 
Emn io EH Trinity 将 clean reads 进行 de novo 
组 装 成 转录 本 , 共 得 到 209 002 条 转录 本 ,平均 长 度 
JJ 687.55 bp, N50 为 982 bp ,序列 长 度 大 于 500 bp 
的 有 87 808 条 , 占 总 序列 数目 的 75. 3496 ( 表 1)。 
对 Trinity 拼装 得 到 的 转录 本 去 元 余 , 共 获得 80 438 
条 unigene ,平均 长 度 为 652.78 bp, N50 为 992 bp, 
其 中 , 28 748 条 序列 长 度 大 于 500 bp , 占 总 序列 数 
目的 35. 74% ; 序列 长 度 在 1 000 bp 以 上 的 有 
13 057 条 , 占 总 数 的 16.23% (R1 和 图 1)。 











表 1 草地 贪 夜 蛾 雄性 成 由 和 $ 龄 幼虫 转录 组 的 组 装 结果 


Table 1 Summary of assembled in the transcriptomes of male adult and the 5th instar larva of Spodoptera frugiperda 








总 数量 平均 长 度 (bp) 
m N50 (bp) N90 (bp) z500 bp 三 1 000 bp TE 
Total number Average length 
转录 本 Transcripts 209 002 982 291 87 808 37 689 687. 55 
独立 基因 Unigenes 80 438 992 264 28 748 13 057 652.78 





2.2 Unigene 功能 注释 

草地 贪 夜 蛾 转录 组 数据 组 装 后 的 unigene 序列 
Hj CDD, KOG, COG, NR, NT, PFAM, Swiss-Prot 
fll TrEMBL 等 多 个 数据 库 比 对 ,共有 46 198 条 
(57. 43% ) unigene 在 至 少 1 个 数据 库 中 获得 功能 注 
释 ,1 713 条 (2. 13% ) unigene 在 所 有 数据 库 中 均 能 
获得 注释 ( 表 2)。 其 中 NR 数据 库 中 注释 到 的 
unigene 数目 最 多 (29 42), unigene 的 
36. 60% ;其 次 为 TrEMBL 数据 库 ,注释 到 的 unigene 
为 29 189 条 , 占 36. 2996 , KEGG 数据 库 中 注释 到 





的 unigene 数目 最 少 (3 408 条 , 4.2496 ) ,涉及 271 
个 代谢 通路 。 尚 有 34 240 条 unigene 没有 得 到 有 效 
注释 , 占 42.57% ( 表 2)。 通 过 与 NR 库 的 比 对 , 获 
得 草地 贪 夜 蛾 unigene 序列 与 近 缘 种 属 的 近似 情况 
并 获得 同 源 序列 的 功能 信息 ,共有 29 393 条 unigene 
获得 注释 。 匹 配 较 多 的 物种 主要 有 脐橙 蜡 蛾 
Amyelois transitella KÆ B. mori , TH Tài PA. Papilio 
xuthus 和 金 凤 蝶 Papilio machaon ,分 别 占 21. 4096 , 
18.75% , 9.3496 和 8. 1796 。 可 见 在 转录 组 水 平 上 
草地 贪 夜 蛾 与 脐橙 蜡 蛾 的 序列 相似 度 最 高 (图 2) 。 
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图 1 草地 贪 夜 峨 雄性 成 由 和 5 龄 幼虫 转录 组 unigene 序列 长 度 分 布 


Fig. 1 


R2 草地 贪 夜 蛾 雄性 成 由 和 S 龄 幼虫 转录 组 
unigene 的 功能 注释 
Table 2 Functional annotation of assembled unigenes 
in the transcriptomes of male adult and the Sth instar 
larva of Spodoptera frugiperda 























F : 基因 数量 
数据 库 百分比 
Number of 
Database Percentage 
genes 
CDD 13 837 17.20 
KOG 13 352 16.60 
NR 29 442 36.60 
NT 28 688 35.66 
PFAM 10 546 13.11 
Swiss-Prot 18 112 22.52 
TrEMBL 29 189 36.29 
GO 23 579 29.31 
KEGG 3 408 4.24 
少 右 一 个 天 数据 诸 中 注释 
至 少 有 一 个 在 数据 库 中 注释 jenos am 
Annotated in at least one database 
"TE 
所 有 数据 库 中 注释 i us 
Annotated in all database 
总 基因 数量 
基因 数量 80 438 100 
Total number of genes 
草地 贪 夜 蛾 的 转录 组 通过 与 NR 数据 库 进行 同 源 比 

















对 和 基因 注释 ,大 量 与 耐 药性 相关 的 基因 家 族 被 鉴 
定 到 ,主要 包括 :310 个 细胞 色素 P450 ( cytochrome 
P450, P450) 基 因 , 72 AA Wt H Bk S- 转 移 酶 
(glutathione S-transferase , GST) 基因 ,161 ^13 BE BS Ij 
(carboxylesterase, COE ) 基因 , 46 个 氨基 上 肽 酶 N 
( aminopeptidase N, APN) 基因 和 91 个 ABC 转运 和 蛋 
白 (ABC transporter) 基因 ( 表 3)。 




















Length distribution of assembled unigenes in the transcriptomes of male adult and the 5th instar larva of Spodoptera frugiperda 


2.3 Unigene 功能 分 类 

将 组 装 得 到 的 草地 贪 夜 蛾 unigene 与 KOG 数据 
库 比 对 ,共有 14 926 条 unigene 获得 注释 , 按 其 功能 
共 分 为 25 类 。 其 中 信号 转 导 机 制 注释 到 的 unigene 
数目 最 多 ,共有 2 000 条 , 占 13.40% ; 仅 适用 于 一 般 
功能 预测 翻译 后 修饰 -蛋白 质 更 新 -分 子 伴侣 .转录 
及 碳水 化 合 物 的 运输 和 代谢 注释 到 的 unigene 占 比 
分 别 为 7.73% , 0. 8996 , 1. 5696 和 4.99%。 仅 有 
48 和 32 条 unigene 分 别 注释 到 核 结构 和 细胞 运动 ， 
分 别 占 比 0. 3296 和 0. 21% (图 3)。 此 外 ,892 条 
unigene 注释 到 未 知 功能 。 由 此 可 见 ,草地 贪 夜 蛾 
unigene 涉及 信号 转 导 功 能 最 多 ,为 以 后 研究 草地 贪 
夜 蛾 代谢 物质 的 调控 等 提供 基因 资源 。 

根据 GO 数据 库 注释 结果 ,共有 205 269 条 草地 
AAIR unigene 被 划分 为 68 个 功能 分 类 (图 4) , 生 
物 学 过 程 注释 到 的 unigene 最 多 ,共有 91 666 个 , 占 
44.66% ,其 中 注释 较 多 的 功能 分 别 为 细胞 过 程 
(14 535 个 ,7.08% ) 代谢 过 程 (11 439 个 , 5.57% )、 
生物 调控 (8 174 个 , 3. 98% ) 和 生物 过 程 调控 (7 258 
个 , 3.53% ) 。 注 释 到 分 子 功能 的 unigene 数目 为 
31 404 个 , 占 15.30% ,其 中 注释 较 多 的 功能 分 别 为 
结合 (13 448 个 , 6. 5596) 、 催 化 活性 (11 108 个 ， 
5.4196 ) 和 转运 体 活性 (1 817 个 , 0.8896) 。82 199 
个 (39. 6196 ) unigene 注释 到 细胞 组 分 ,其 中 注释 较 
多 的 功能 分 别 为 细胞 (14 977 个 , 7. 2996) 、 细 胞 部 
分 (14 934 个 , 7.27% ) .细胞 器 (11 569 个 , 5. 63% ) 、 
膜 (9 401 个 , 4.58% ) 和 腊 部 分 (7 534 个 , 3.67% ) 。 
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图 2 ”草地 贪 夜 峨 雄性 成 了 由 和 5 龄 幼虫 转录 组 unigene 在 NR 数据 库 的 同 源 物 种 分 布 
Fig. 2 Species distribution of the top BLAST hits against the NR database for the assembled unigenes 





in the transcriptomes of male adult and the 5th instar larva of Spodoptera frugiperda 


表 3 草地 贪 夜 蛾 雄性 成 由 和 S5 龄 幼虫 转录 组 中 与 解毒 相关 基因 的 统计 
Table 3 Statistics of unigenes associated with detoxification in the transcriptomes 


of male adult and the 5th instar larva of Spodoptera frugiperda 


























抗 药性 基因 家 族 总 基因 数量 数据 库 Database 
Insecticide tolerance Total number 
gene families of genes NR PFAM Swiss-Prot 
P450 310 310 193 263 
GST 72 67 25 58 
COE 161 161 142 141 
APN 46 45 28 37 
ABC 91 84 33 66 


x10 





A: RNA 加 工 和 修饰 RNA processing and modification 

B: 染色 质 结构 和 动力 学 Chromatin structure and dynamics 

C: 能 源 生产 和 转化 Energy production and conversion 

D: 细胞 周期 控制 、 细 胞 分 裂 、 染 色 体 分 区 Cell cycle control, cell division, chromosome partitioning 
E: 氨基 酸 转 运 和 代谢 Amino acid transport and metabolisn 

F: 核 背 酸 转运 和 代谢 Nucleotide transport and metabolism 

G: 碳水 化 合 物 的 运输 和 代谢 Carbohydrate transport and metabolism 

H: 辅酶 转运 和 代谢 Coenzyme transport and metabolism 

1: 脂 质 运输 和 代谢 Lipid transport and metabolism. 

上 上 翻译、 核糖 体 结构 和 生物 发 生 Translation, ribosomal structure and biogenesis 
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5 1.254 K: 转录 Transcription 
mz L: 复制 、 重 组 和 复制 Replication, recombination and repair 
B E 1.00 M: 细胞 壁 / 膜 / 包 膜 的 生物 发 生 Cell wall/membrane/envelope biogenesis 
E E : N: 细胞 运动 Cell motility 
o O: 翻译 后 修饰 、 蛋 白质 更 新 、 分 子 伴侣 Posttranslational modification, protein turnover, chaperones 


P: 无 机 离子 的 运输 和 代谢 Norganic ion transport and metabolism 

Qi: 次 生 代谢 产物 的 合成 、 转 运 和 分 解 代谢 Secondary metabolites biosynthesis, transport and catabolism 
R: 仅 适 用 于 一 般 功 能 预测 General function prediction only 

S: 仅 能 未 知 Function unknown 

T: 信号 转 导 机 制 Signal transduction mechanisms 

U: 细胞 内 运输 、 分 泌 和 圳 泡 运输 Intracellular trafficking, secretion, and vesicular transport 

V: 防御 机 制 Defense mechanisms 

W: 细胞 外 结构 Extracellular structures 

Y: 核 结构 Nuclear structure 

Z: 细胞 骨架 Cytoskeleton 
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图 3 草地 贪 夜 蛾 雄性 成 由 和 5 龄 幼虫 转录 组 unigene 的 KOG 功能 分 类 


Fig. 3 KOG function classification of unigenes in the transcriptomes of male adult and 





the 5th instar larva of Spodoptera frugiperda 
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图 4 草地 贪 夜 蛾 雄性 成 由 和 5 龄 幼虫 转录 组 unigene 的 GO 功能 分 类 


Fig. 4 GO classification of unigenes in the transcriptomes of male adult and the 5th instar larva of Spodoptera frugiperda 





根据 KEGG 数据 库 注 释 结果 , 主要 分 为 5 大 类 
功能 ,包括 代谢 (1 815 个 ，29. 25% )、 有 机 系统 
(1 496 个 , 24. 11% ) 遗传 信息 处 理 (1 118 个 ， 
18. 0196 ) ` 细 胞 过 程 (917 个 ,14.78% ) 和 环境 信息 
处 理 (860 个 , 13.8696) (图 5)。 在 代谢 中 ,与 氨基 
酸 代 谢 和 脂 质 代 谢 相 关 的 unigene 分 别 为 275 个 
(4.43% ) 和 249 个 (4.01% ) 。 在 有 机 系统 中 ,涉及 
环境 适应 的 unigene 有 68 个 , 占 1.09% 。 在 遗传 信 
息 处 理 中 注释 到 unigene 最 多 的 代谢 过 程 为 翻译 
(433 个 , 6.97% ) , 而 在 细胞 过 程 中 ,与 运输 和 分 解 
代谢 相关 的 unigene 最 多 (379 个 , 6. 11% ) 。 信 和 号 
转 导 在 植物 发 育 过 程 以 及 对 外 界 刺 激 的 响应 过 程 起 
重要 作用 ,本 研究 共 检 测 到 739 条 unigene 涉及 信和 号 
转 导 ,为 下 一 步 鉴定 调控 草地 贪 夜 峨 生长 发 育 的 信 
号 因子 提供 研究 基础 。 

2.4 草地 贪 夜 蛾 雄性 成 虫 与 $ 龄 幼虫 差异 表达 基 
因 分 析 

雄性 成 贝 有 69 204 个 基因 表达 ,5 龄 幼虫 有 
65 634 个 基因 表达 。 雄 性 成 虫 和 5 龄 幼虫 共同 表达 
基因 有 54 990 个 ,分 别 占 雄性 成 虫 转录 组 的 
79. 46% 和 5 龄 幼虫 转录 组 的 83.78% 。 进 一 步 鉴 
定 得 到 雄性 成 虫 阶段 特异 表达 的 基因 有 14 214 个 ， 
5 龄 幼虫 阶段 特异 表达 的 基因 有 10 644 个 (图 6: 
A)。 以 5 龄 幼虫 为 对 照 组 ,分 析 草 地 贪 夜 峨 2 个 不 


Av 





同 生长 发 育 阶 段 基因 表达 的 差异 ,在 雄性 成 虫 与 5 
AFA n cc E de 4E QE 15 561 个 差异 表达 基因 
( differentially expressed genes, DEGs) ,其 中 9 162 个 
上 调 表达 ,6 399 个 下 调 表 达 ( 图 6: B). GO 功能 富 
集 分 析 发 现 ,在 上 调 表达 基因 中 鉴定 到 3 个 基因 涉 
及 信息 素 结合 (pheromone binding) (GO:0005550) , 
1 个 基因 涉及 信息 素 代 谢 过 程 (pheromone metabolic 
process) (GO:0042810) 以 及 大 量 涉 及 脂肪 酸 合 成 和 
延长 等 功能 的 基因 ,这 些 在 雄性 成 虫 中 上 调 表达 的 
基因 为 下 一 步 鉴定 草地 贪 夜 蛾 性 信息 素 合 成 提供 了 
可 能 。KEGG 功能 富 集 发 现 雄 性 成 虫 上 调 表 达 的 基 
中 涉及 信号 转 导 的 通路 显著 富 集 ,例如 : 胰 高 血糖 
素 信 号 通路 ，AMPK 信号 通路 ( AMPK signaling 
cGMP-PKG 信号 通路 ( cGMP-PKG 
signaling pathway ) 和 钙 信 号 通路 (calcium signaling 
pathway) (图 7: A) 。 在 雄性 成 虫 下 调 表 达 基 因 中 , 
我 们 发 现 涉 及 解毒 相关 的 通路 ,例如 :细胞 色素 
P450 对 异 生物 的 代谢 (metabolism of xenobiotics by 
cytochrome P450) 药物 代谢 -细胞 色素 P450 (drug 
metabolism-cytochrome P450 ) 和 药物 代谢 -其 他 酶 
(drug metabolism-other enzymes) 显著 富 集 ( 图 7: B), 
2.5 草地 贪 夜 蛾 嗅觉 相关 基因 分 析 

在 本 研究 中 ,通过 同 源 比 对 , 共 鉴 定 到 424 ^ Dt 
觉 相 关 的 基因 ,包括 38 个 化 学 感应 蛋白 
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成 虫 和 5 龄 幼虫 转录 组 unigene 的 KEGG 功能 分 类 


Fig. 5 KEGG classification of unigenes in the transcriptomes of male adult and the 5th instar larva of Spodoptera frugiperda 
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草地 贪 夜 蛾 雄性 成 由 和 5 龄 幼虫 差异 表达 基因 的 比较 
Fig. 6 Comparison of differentially expressed genes ( DEGs) in male adult and the 5th instar larva of Spodoptera frugiperda 


A; 雄性 成 虫 与 5 龄 幼虫 基因 表达 差异 的 维 恩 图 Venn diagram showing the DEGs between male adult and the 5th instar larva; B; 雄性 成 里 与 龄 幼 


虫 差异 表达 基因 的 散 点 图 Scatter plot showing the DEGs between male adult and the 5th instar larva. TPM; 每 


million. 
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图 7 草地 贪 夜 蛾 雄性 成 虫 与 5 
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Fig. 7 KEGG pathway enrichment of the differentially expressed genes ( DEGs) in male adult 





A: 雄性 成 虫 中 比 5 


and the 5th instar larva of Spodoptera frugiperda 











adult and the 5th instar larva; B; ffl 











regulated DEGs between male adult and the 5th instar larva. 


(chemosensory protein) 基因 、48 个 气味 结合 蛋白 
(odorant binding protein ) 基因 、294 个 气味 受 
(odorant receptor ) 基因 、14 个 嗅觉 受 体 (olfactory 


receptor) 基因 和 30 个 味 党 受 


基因 ( 表 4)。 进 一 步 分 析 基 因 表 达 的 差异 发 现 ,在 


体 


觉 受 体 ( gustatory receptor ) 


雄性 成 由 和 5 龄 





龄 幼虫 中 ,涉及 嗅觉 受 
差异 表达 ,而 其 他 几 类 嗅 党 相关 基因 中 ,在 雄性 成 虫 
中 上 调 表达 的 数量 明显 多 余 5 龄 幼虫 ,这 可 能 与 雄 
性 成 虫 羽 化 后 的 各 种 生理 活动 有 关 , 因 此 这 些 基因 
可 能 作为 草地 贪 夜 蛾 防 控 RNA 干扰 的 候选 靶 标 。 


龄 幼虫 上 调 表 达 基 因 的 KEGG 代谢 通路 富 集 分 析 KEGG pathway enrichment analysis of the up-regulated DEGs between male 
成 虫 中 比 5 龄 幼虫 下 调 表 达 基 因 的 KEGG 代谢 通路 富 集 分 析 KEGG pathway enrichment analysis of the down- 


体 的 基因 没有 
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表 4 草地 贪 夜 蛾 雄性 成 由 和 $ 龄 幼虫 嗅觉 相关 基因 的 表达 分 析 


Table 4 Expression analysis of olfaction-related genes in male adult and the 5th instar larva of Spodoptera frugiperda 























总 注释 数 差异 表达 基因 
嗅觉 相关 基因 Total number of Differentially expressed genes (DEGs) 
Olfaction-related genes annotated 总 数 Es 下 调 

unigenes Total Up-regulated Down-regulated 
化 学 感应 蛋白 基因 Chemosensory protein genes 38 19 13 6 
气味 结合 蛋白 基因 Odorant binding protein genes 48 32 24 8 
气味 受 体 基因 Odorant receptor genes 294 55 38 17 
嗅觉 受 体 基 因 Olfactory receptor genes 14 0 0 0 
味觉 受 体 基因 Gustatory receptor genes 30 3 3 0 

个 基因 编码 CST,161 个 基因 编码 COE ,46 个 基因 编 
3 讨论 码 APN 和 91 个 基因 编码 ABC 转运 蛋白 ( 表 3)。 




















近年 来 , 随 着 高 通 量 测序 技术 的 发 展 ,凭借 其 产 
出 数据 量 大 、 成 本 低 和 产 出 速度 快 等 优势 ,在 新 基因 
发 据 、 功 能 基因 鉴定 和 分 子 标记 开发 上 的 应 用 越 来 
越 广 泛 。 在 非 模式 昆虫 的 研究 中 ,例如 其 他 鲜 翅 目 
Eb rk d £X $E Spodoptera litura, 小菜 蛾 Plutella 
xylostella 等 ,高 通 量 测序 技术 已 经 成 为 研究 其 功能 
基因 的 挖 气 、 遗 传 进化 等 方面 的 有 效 手段 (Vogel et 
al., 2011; Gu et al., 2013; You et al., 2013 ) 。 本 研 
究 基 于 Illumina HiSeq2500 的 高 通 量 测序 平台 对 草 
地 贪 夜 蛾 雄性 成 虫 和 5 龄 幼虫 进行 de novo 转录 组 
测序 , 共 获 得 80 438 条 unigene ,平均 长 度 为 652. 78 
bp, N50 为 992 bp, 

基因 功能 注释 结果 显示 , 共有 46 198 条 
(57. 4396 ) unigene 在 至 少 一 个 数据 库 中 获得 功能 注 
释 ,1 713 条 (2. 13% ) unigene 在 所 有 数据 库 中 均 能 
获得 注释 ,但 仍 有 约 42.75% 的 unigene 没有 获得 相 
关注 释 信 息 , 可 能 是 由 于 部 分 序列 组 装 长 度 过 短 , 缺 
少 保守 的 核心 序列 无 法 获得 准确 功能 注释 ,也 可 能 
是 非 编码 序列 或 者 是 功能 未 知 的 新 基因 ( Hou. et al., 
2011)。 在 KOG 数据 库 和 KEGG 数据 库 中 分 别 鉴定 
到 2 000 和 739 条 unigene 涉及 信号 转 导 功 能 ,为 下 
一 步 解 析 草 地 贪 夜 蛾 响应 外 界 刺激 提供 研究 基础 。 

在 杀 虫 剂 长 期 使 用 过 程 中 ,昆虫 容易 产生 不 同 
程度 的 耐 药性 。 高 通 量 测序 技术 能 够 从 基因 层面 解 
析 昆 虫 的 耐 药 性 和 适应 性 机 理 , 为 研发 农药 的 杀 虫 
选择 靶 点 提供 参考 。Cheng 等 (2017 ) 比较 基因 组 分 
析 发 现 斜 纹 夜 蛾 中 涉及 解毒 相关 的 基因 家 族 
(P450, GST, COE, APN 和 ABC) 急剧 扩 张 。 在 本 
研究 中 ,我 们 在 草地 贪 夜 蛾 转录 组 中 鉴定 到 大 量 涉 
及 解毒 相关 的 基因 ,包括 310 个 基因 编码 P450 ,72 









































































































































昆虫 通过 高 度 灵敏 和 复杂 的 嗅觉 系统 感知 环境 
中 的 化 学 信息 来 完成 一 系列 生理 活动 ,因此 对 昆虫 
的 嗅觉 系统 进行 干扰 是 防治 虫害 的 有 效 途 径 
(Plettner et al., 2002; Will and van Bel, 2006), Eh 
虫 嗅觉 过 程 需 要 多 种 与 嗅觉 相关 和 蛋白 及 受 体 参 与 ， 
153] ^ (2016 ) 利用 高 通 量 技术 在 荔枝 蒂 星 虫 
Conopomorpha sinensis 转录 组 中 挖掘 到 100 个 涉及 
嗅觉 相关 的 基因 。 本 研究 中 我 们 共 鉴 定 到 424 个 基 
因 涉及 嗅觉 相关 基因 ,包括 38 个 化 学 感应 蛋白 基 
48 个 气味 结合 蛋白 基因 、294 个 气味 受 体 基 因 、 
14 个 嗅觉 受 体 基因 和 30 个 味觉 受 体 基因 ( 表 4) ， 
并 且 在 雄性 成 由 和 5 龄 幼虫 之 间 的 表达 存在 差异 ， 
表明 雄性 成 虫 与 5 龄 幼虫 在 嗅觉 感知 方面 存在 差 
异 ,可 能 与 不 同 发 育 阶 段 的 生理 活动 差异 有 关 LONE 
发 针对 草地 贪 夜 蛾 不 同 发 育 阶段 特异 性 农药 提供 新 
RE. 
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